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^ I give a physical explanation of what shear viscosity is, and what physics de- 

, termines its value. Then I explain why determining the shear viscosity of the 

' Quark-Gluon Plasma is interesting. I outline the leading-order calculation of the 

QQ ' QGP shear viscosity (and baryon number diffusion constant), explaining why the 

^ quite complicated physics of parton splitting and Landau- Pomeranchuk-Migdal in- 

, terference effects are required for its calculation. Then I briefly explore the range 

' of applicability, emphasizing the importance of plasma instabilities. 

' , 1. Introduction 

Ph. 

1^ ■ In this talk I will discuss work with my collaborators, Peter Arnold and 

Larry Yaffe, on transport coefficients in QCD, particularly shear viscosity. 

Interest in transport coefficients in the Quark-Gluon Plasma (QGP) 
started in the 1980's, and several rough estimates were made at that time 



^ I ^. The first really mature work was by Baym et al. ^ in 1990. They 

showed that the dominant physics involved is small angle scattering, and 
they demonstrated how to perform a calculation accurate to leading order 
in the logarithm of the strong coupling, \n{l/gs). We started work around 
1999 under the mistaken impression that not much was required to improve 
their work to leading order in the coupling gg. After discovering that even 
the leading-log calculation was wrong in the literature, we decided to write 
two papers, rather than one-a quick paper ^ computing the leading log 
correctly, and a longer one presenting a leading order calculation. The 
quick paper took 8 months and the "second" paper turned out to be six 
papers 4,5,6,7,8,9^ spanning the next two and a half years. 

My goal is to explain how the physics of transport coefficients turns out 
to be much richer than we had anticipated. In Section 2, I will explain 
what shear viscosity is and what physics is involved; and I will briefly 
explain why it is an interesting property of the QGP. Section 3 will outline 
our calculation, emphasizing the roles played by soft scatterings, identity 
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changing scatterings, and collinear processes. Then. Section 4 will briefly 
discusses how limited the range of the calculation is; not only does it depend 
on weak coupling, but on a near-equilibrium assumption outside of which 
the physics is completely different. 

2. The physics of shear viscosity 

The velocity of an ideal fluid evolve according to Euler's equations, 

dv _ 1_ , , 

— = -vVv - -Vp, P = p{p). (1) 

Furthermore, the velocity at a boundary equals the velocity of that bound- 
ary (no slipping) , and the stress tensor is determined by the pressure alone, 
that is, in the v = frame it is given by 

T^'^ = {p+ p)u''u'' + prj'"'' , (2) 

(with 7]^^" — Diag[— , +, +, +]). In fact, Euler's equations follow simply 
from the assumption that the stress tensor will satisfy this form, from an 
equation of state, and from stress-energy conservation, df^T'^'^ = 0. 

Local equilibrium ensures that these conditions are met. What enforces 
local equilibrium is interactions between degrees of freedom; free theories 
never equilibrate. Therefore, the ideal description is a better and better 
approximation as a theory is more and more strongly coupled; in this case 
ideality is the opposite of free behavior. 

To see this, consider a fluid undergoing shear flow. Shear flow means 
that, in the local rest frame, the traceless part of diVj is nonzero. The 
simplest example of shear flow is a fluid trapped between two plates, one 
flxed and one moving laterally with velocity v, as in Fig. 1; under ideal 

/^ / ^ / 7^ /^/ / t:^" 

^ ^ ^ ^ Flow varies 

with height 



Fixed bottom 

/////////// / bo^^dary 
Figure 1. a system undergoing shear flow 



Eulerian flow, the fluid at the top is moving with the upper boundary, with 
Txz = 0; therefore there is no lateral force on the top boundary. 
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In a real fluid, there will be a "drag" force on the top plate; the shear 
viscosity is defined by the relation that the force per unit area is, 

^ = Tlrz = , implying = pSij +r] (^diVj+djVi~SijdkVk^ . (3) 

Shear occurs because particles can fly freely. The particles near the top 
of the fluid region arc primarily moving to the right; those near the bottom 
are moving right or left with equal likelihood. If the particles fly freely, the 
upward moving particles near the bottom, and the downward moving parti- 
cles from the top, will meet in the middle, as shown in Fig. 2. However, the 




Figure 2. Free flow leads to an asymmetrical momentum distribution. 



particles coming down from the top tend to be moving to the right, while 
particles coming up have no preferred direction; so the momentum distri- 
bution in the center (and everywhere, actually) tends to become skewed. 

Such a skewed distribution has nonvanishing T^z and there is net transverse 
momentum flow and a net drag on the top surface. Scatterings restore the 
distribution towards equilibrium, evening out the momentum distribution. 
The fewer scatterings, the more distorted the momentum distribution will 
be. Therefore, the shear viscosity will be smaller if the scattering rate is 
larger; 77 oc l/F with F the mean scattering rate. ^ The problem of deter- 

Familiar viscous fluids like honey are viscous because the molecules get tangled in each 
other, and are not a good analogy with any system which behaves approximately like a 
collection of freely flying objects. A better example is the air; as the pressure is lowered, 
there are fewer particles available to carry momentum, Tij oc p, so one would naively 

expect rj oc p. However, at lower pressure, collisions are rarer, T oc p; and so the shear 
viscosity of air is actually independent of air pressure. 
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mining the shear viscosity will be the problem of determining how collision 
processes re-arrange the momentum distribution of particles from skewed 
towards isotropic. 

Why is the shear viscosity of the QGP interesting? One of the main 

focuses of this conference is the understanding of nonequilibrium dynamics; 
QCD is the most interesting theory to study this and viscosity and other 
transport coefficients (baryon number diffusion and electrical conductivity, 
which have very similar physics) are theoretically clean and well defined 
nonequilibrium dynamical quantities. That makes them ideal objects for 
"cutting our teeth" on calculational strategies for nonequilibrium field the- 
ories. 

The shear viscosity of QCD is also directly interesting because it gives 
a handle on the quality of the hydrodynamical description of heavy ion 
collisions; this has been considered at some length recently by Teaney 
who found the corrections to flow parameters like V2 due to viscosity at 
RHIC. Viscosity and particularly diffusion coefficients are also important 
in electroweak baryogenesis. 



3. The shecir viscosity of the QGP 

Shear viscosity describes momentum transfer, and most of the momentum 
of a weakly coupled, near-equilibrium QGP is in quasiparticle excitations 
with momentum p ^ 3T. These should be adequately described by kinetic 
theory, as I implicitly assumed in the last section's discussion. In fact, the 
validity of a kinetic description has essentially been derived diagrammati- 
cally, by Jeon for At/)'* theory and by others for gauge theories i2,i3^ 

The physics of viscosity is the physics of direction changing scattering, 
and this is dominated by 2 2 f-channel gluon exchange ^. In terms of the 
spatial momentum q exchanged in the scattering process, such scatterings 
have a cross section of form, 

da(x J , (4) 

q ^4 

which is quadratically small q divergent — the familiar Coulomb divergence. 
However, what is important for shear viscosity is how quickly the parti- 
cles change direction. Momentum changes from different scatterings add 
incoherently; so 



(5) 



February 2, 2008 0:46 Proceedings Trim Size: 9in x 6in 



moore 



is the relevant quantity and enough coUisions must occur for it to add up 
to for a particle to change direction. This integral is log IR divergent. 
There are two key lessons. 

(1) Small q collisions dominate, but only by a logarithm. 

(2) Large momentum particles take the longest time to change direc- 
tion; since (Sp)'^ ^ is required, the time scale for direction 
change goes as r oc Therefore high energy particles contribute 
disproportionally to rj. 

Baym et al? provided the first quantitative treatment which accounted 
for both of these points. To find the shear flow, they wrote the Boltzmann 
equation for the nonequilibrium distribution function f{x,p), 

g/' ^^ =-v- V./b, X, t] - C[f] [p, X, t] , (6) 

and linearized it. That is, they took f = f^ + Sf and expanded the collision 
integral C to flrst order in 6f. To this order, one can write 

Sf = /o(l ± fo)diVjipiPj - ^)x(p) . (7) 

The collision integral is nonlocal in momentum space, so we have an integral 
equation to solve, which in general is nasty. They showed how to approach 
this variationally. They also showed that the physics of dynamical screen- 
ing by the plasma, which is contained within the Hard Thermal Loops 
(HTL's) of Braaten and Pisarski is suflacient to render the Coulomb 
divergence finite; the would be logarithmic divergence becomes a large log- 
arithm ln{l/g). Further, at leading order in this logarithm C is semilocal 
in momentum, which allowed Baym at al. to solve for the shear viscosity 
of QCD at leading logarithmic order. 

In our first paper ^, we showed that, even at leading logarithmic order, 
one must also include Compton type processes, 






which are logarithmically IR divergent. The role of these processes is not 
to change particle direction, but to change particle identity; a quark on 
scattering becomes a gluon and vice versa. Because of their different group 
Casimir, gluons equilibrate more quickly than quarks; so the conversion of 
hard to equilibrate quarks into easy to equilibrate gluons accelerates ther- 
malization. This is needed at leading-log order, though it is only numeri- 
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cally important for electrical conductivity in QED. With these processes, a 
leading- log treatment is possible ^, yielding, for instance, 



3n'^ + 32Nspec e2 1nl/e' 



(8) 



_ 2^3^C^(3)7r-3 1 

The result for rj is more cumbersome and can be found in our paper ^ . The 
above results are only accurate to 0.5% because they were made with the 
Ansatz ^(p) cx p^. More accurate results are presented in the paper. 

What is required to convert these to a leading-order computation? One 
only needs to treat the collision integral more carefully, both in the q ^ T 
momentum region where the collision integral is nonlocal, and in the q ~ 
Tod ~ 9sT region where screening effects must be dealt with. For heavy 
quark diffusion, this is sufScient; Derek Teaney and I recently derived a 
complete leading order result; 

(10) 



TVr \ f, 2T 1 C'(2)\ Nf In 2 

^c+^ In + 7:-7i3+ ' ' 





CfQlT ' 2 ; \'"mD '2 ■ C(2) 

However, viscosity in the QGP and light quark diffusion turn out to in- 
volve more physics, and therefore more work; one must also include inelastic 
processes with collinear external states. 



Particles in the plasma undergo the left process at an 0{asT) rate, not a^T, 
because of the Coulomb divergence already discussed. The right process is 
slower by a factor of as; this is true even when the radiated gluon has 
momentum 0{T), provided it is collinear with the particle which emits it. 
The rate for this process is therefore as large as the large angle scattering 
rate. 

We saw previously that the particles furthest from equilibrium were 
those with the most energy, 6f oc f\p\^. Collinear processes like the one 
just presented tend to split these high energy particles into lower energy 
ones, which are more easily bent; therefore they lead to equilibration and 
lower 77, even though they do not change particle direction. Numerically, 
we find that they reduce rj by about 10% ^. 

But life is even more complicated; one cannot ignore interference be- 
tween emission amplitudes of form 
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The effect of sucli interference is referred to as the Landau-Pomeranchuk- 
Migdal (LPM) effect , and it suppresses the efficiency of collinear emission. 
To see the physical origin of the suppression, consider a particle undergoing 
multiple scatterings, and remember that it possesses a finite transverse wave 
packet. That wave packet continues to physically overlap the wave packet 
of the emitted radiation for some time, called the formation time; if there 
is a second scattering in that time, the two radiation fields overlap and 
interfere; 



Second scattering - 



Approaching quark 




Second emitted 
photon 



First scattering 

The radiations sum in amplitude; on average the interference is destructive. 

Dealing with this complication takes a lot of work, and has been the 
topic of several studies ^'^•^•'^ . The upshot is that the rate at which a particle 
of momentum p undergoes a soft scattering which induces a bremsstrahlung 
of a gluon of energy fc, is 



dTjp, k ) _ Cs9l 
dkdt 



1 



1 



1 ± e-fe/r 1 ± e-(v-k)/T 



qq 
99 



:2h-ReF{h,p,k) , 



> X 



(11) 



where F is the solution of the following integral equation: 



2h = iSEF{h) + g 



F{h)-F{h-kq^)] 



Ca 
' 2 



5E 



2pk{p-k) 2k 



with C(qj 



[F{h) - F{h+pq^) + Fih) - F{h-{p-k) g^)]}, 

(12) 

(2iVc+^f). Determining dV/dkdt 



'^p-k _ % 
2{p-k) 2p 

™2 _ a^T^ 
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therefore requires the solution of an integral equation; and this for a collision 
term which, to solve the Boltzmann equation, must appear as a coefficient 
in an integral equation which is to be solved. However, the numerics turn 
out not to be very demanding. 

These ingredients have been combined into a kinetic theory for the QGP 
which has been linearized and solved to determine the transport coeffi- 
cients at leading order in the strong coupling ^. Results for shear viscosity 
are presented in Fig. 3. At ckg = 1/3, r/ is quite small, and hydrodynamics 
should be nearly ideal; but the perturbative computation is breaking down 
at this point and the results cannot be trusted (factor of 3 uncertainty?). 




Coupling strength a. 



Figure 3. Shear viscosity, scaled by Og ^ (left) and by enthalpy density (right). 



4. Range of applicability of QGP transport coefficients 

Two key assumptions underlie the above calculation. The first is that the 
coupling is actually weak, at the scale relevant to the physics in question. 
For electroweak phase transition applications this is probably a good ap- 
proximation; for RHIC it is dubious, though most of the theory community 
is guilty of this assumption. I will not discuss this more here. 

The other assumption underpinning our treatment is, that the system 
is near equilibrium. One can check self-consistency by comparing the time 
scale of relaxation with the age of the plasma. The ratio 'ri/{p + p) is a time 
or length scale, called the sound attenuation length; it differs by a geometric 
factor ~ 1/5 from the relaxation time of the plasma. At as — 1/3, Fig. 3 
shows that this is ^ 2/3T, so r ~ S/T, which is a few fermis at RHIC. 
Therefore we expect hydrodynamic treatments to be valid after ^ 3 fermi 
at RHIC. One might question a hydro description for times earlier than 
that. 
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How different could the early time physics be? The answer turns out 
to be, pretty different. For systems which are far from equilibrium but ap- 
proximately isotropic, kinetic theory should be applicable under fairly weak 
conditions ^. But for systems with strong momentum space anisotropy, this 
turns out not to be the case. A treatment involving soft gauge fields be- 
comes essential, because these soft gauge fields show exponential growth 
due to plasma instabilities It has recently been noted that this may be 
important for QGP physics i9'20,2i^ ^^^^ jg ^^iq topic of two other talks 
at this conference 22,23 r^j^g short story is that, for plasmas which are far 
from isotropic, soft gauge field modes can grow at a rate suppressed by only 
a single power of Qs, and may randomize the directions of other particles 
on this time scale. 



5. Conclusions 

Transport coefficients in gauge theories can be computed within kinetic 
theory, but they involve rich physics; hard thermal loops, multiple scatter- 
ing processes, coUinear splitting. The kinetic description is only possible 
after resumming certain multiple scattering emission processes to account 
for LPM interference. It has taken nearly 20 years to go from the first at- 
tempts at a quantitative analysis to one which is complete at leading order 
in the coupling g^. 

The calculation of transport coefficients has been directly useful to the 
heavy ion community ^'^j but it has also been a good object lesson in non- 
equilibrium field theory in QCD. Most importantly, we have learned just 
how narrow the domain of validity of kinetic theory is, and how much 
more complicated the physics of highly anisotropic plasmas may be. The 
problem of understanding plasma instabilities in QCD is only beginning to 
be addressed and should be the subject of much future work. 
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